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Spectral Functions, Maximum Entropy Method and Unconventional 
Methods in Lattice Field Theory 

Chris Allton, Danielle Blythe and Jonathan Clowser 

Department of Physics, University of Wales Swansea, Singleton Park, Swansea SA2 8PP, United 
Kingdom 

We present two unconventional methods of extracting information from hadronic 2-point functions produced 
by Monte Carlo simulations. The first is an extension of earlier work by Leinweber which combines a QCD Sum 
Rule approach with lattice data. The second uses the Maximum Entropy Method to invert the 2-point data to 
obtain estimates of the spectral function. The first approach is applied to QCD data, and the second method is 
applied to the Nambu-Jona-Lasinio model in (2+l)D. Both methods promise to augment the current approach 
where physical quantities are extracted by fitting to pure exponentials. 



1. Introduction 

The conventional approach of studying lattice 
hadronic correlation 2-point functions, G2(t), is 
to fit the large time, t > t m i n , interval of G2(t) to 
an exponential, Ze~ Mt , where M is the mass of 
the asymptotic state, and Z is the overlap of that 
state with the interpolating operator used. This 
method can be extended by including further ex- 
ponentials, Zie~ Mit , which correspond to the low- 
est excited states of the channel in question. How- 
ever, fundamentally, the small time part of the 
correlation function, t < t m i n , is always thrown 
away. This paper presents two approaches which 
are able to include the f « region of G2(t) in 
the analysis. 

The first follows the approach of Leinweber |Q 
where a QCD Sum Rule technique is used to 
model the short time behaviour of G2(t). This 
is detailed in the next section where new results 
for non-degenerate mesons are presented, and in 
Section 3 where these calculations are extended to 
include lattice artefacts coming from the Wilson 
action. 

The second approach uses the "Maximum En- 
tropy Method" of numerical analysis to invert 
G 2 (t) into its spectral function representation. 
Section 4 describes this method and its ap- 
plication to the Nambu-Jona-Lasinio model in 
(2+l)D. 



2. QCD Sum Rule approach and non- 
degenerate mesons 

In a method was presented which combines 
two of the main non-perturbative approaches to 
strongly coupled field theory: the lattice and 
QCD Sum Rules. This section presents an 
overview of this approach and its extension to the 
case of non-degenerate mesons with non-zero mo- 
mentum. 

The basic quantity considered is the short dis- 
tance expansion (OPE) of the quark propagator, 

Sq(x), 
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where m is the quark mass. 

The 2-point correlation function at zero mo- 
mentum, C?2(i), can be written as, 



G r 2 (t) = ^<0|T{Jr(x,i)Jr(0,0)}|0>, 

X 

or, Wick contracting, as 



G£(t) = / d 3 xTr{S q (x)TS q (- 



Or}. 



(2) 



(3) 



r is a spinor matrix which defines the quantum 
numbers of the interpolating operator, Jr, and 
thereby of the channel considered. In Eq.(||) we 
have restricted our attention to mesons (although 
in |l| baryons were considered). 



2 



The QCD Sum Rules approach (and the 
method of B) substitutes the short distance 
expansion S® PE into the expression for G2 in 
Eq.(|). This defines G$ PE . 

We now introduce the spectral density, p(s), 
which is defined via the Laplace transform of 
G 2 (t), 



(4) 



,OPE 



corresponding to 



G 3 (t) = / p(s)e- st ds. 
Jo 

The spectral density, p(s 
G° PE , is obviously valid for large values of s 
(i.e. short distances). We cannot expect p(s) OPE 
to reproduce the non-perturbative behaviour of 
QCD which generates bound states. The QCD 
Sum Rules "Continuum Model" allows for this by 
introducing a "continuum threshold", sq, above 
which p(s) is assumed valid. The non- 
perturbative, confining features of the system are 
modelled with a delta function at s — M < s$. 
The full spectral function is, 



p(s)=Z5(s-M) + Z6(s-s )p(s) 



OPE 



(5) 



Z is the normalisation of the ground state whose 
mass is M . and £ a parameter which takes ac- 
count of the different normalisations of the lattice 
versus continuum states (see [Q). 

For completeness, we write here the two-point 
function that Eq.(||) defines: 



G 2 (t) = Ze- Mt + U p(s) 



\OPE n -st 



ds. 



(6) 



The above method was first employed in the 
study of lattice data in jl| where the analytic ex- 
pressions for Eq.(||) were derived in the case of 
zero-momentum baryonic channels. In Q), this 
work was extended to zero-momentum mesons 
and the continuum limit of the Ansatz, Eq.(|^), 
was studied. In ||, we extend this approach 
still further by studying non-degenerate mesons 
at both momentum p = and (1, 1,1). We 

show here the expression for G2 PE in the case of 
the (spatial) vector channel at zero momentum. 
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Parameter 


Value 





6.0 


N f 


(i.e. Quenched) 


Volume 


16 3 x 48 


csw 


(i.e. pure Wilson Action) 



8ttV 
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Table 1 

The lattice parameters used. 



where mi 2 are the two valence quark masses of 
the meson, and the superscript cord refers to the 
fact that the continuum expression for S q , Eq.(|l]), 
was used in the derivation of Eq. ((?]) . 

In this work the Monte Carlo data for G2 (t) was 
generated using UKQCD computing resources. 
Table 1 lists the lattice parameters used in this 
simulation. 

In Figure [l] we show the effective mass for the 
vector channel at K1.2 = 0.1385. In this plot the 
QCD Continuum Model fitting function Eq.(0) is 
shown together with the result of a 2-exponential 
fit. The fitting range for both Eq.(^) and the 2- 
cxponential fit was t = 2 — 20. 

As can be seen from this plot, the QCD Con- 
tinuum Model fit, Eq. ([?]), clearly gives a better 
fit. This can be seen quantitatively in Table || 
where the y 2 values for both fits are shown to- 
gether with the estimates of the Z and M values 
for a wide range of k values. The results of a 1- 
exponential fit (using t = 12 — 20) are shown as 
a control. 

3. Analytic lattice results for non- 
degenerate mesons 

In Section 2, an analytic expression for the 2- 
point function G2 was derived by using the con- 
tinuum, short-distance expansion for the quark 
propagator, S q . In this section we extend this 
work still further using the lattice expression 
for (the co-ordinate space) quark propagator, 
S q (x) LAT , recently derived by Paladini for Wil- 
son quarks Q. 

While the momentum space expression for lat- 
tice propagators has a simple, closed form, the 
co-ordinate space representation does not. Be- 
fore the work of jl, S q (x) LAT was known in two 
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K = 


0.1550 




1-exp 


6.7( 5)xl0" 3 


.426 ( 4) 




2-exp 


8.5( 2)xl(T 3 


.440 ( 2) 


160/15 


CM 


5.9( 2)xl0~ 3 


.419( 3) 


30/15 


0(a) 


5.9( 2)xl0- 3 


.419( 3) 


30/15 




K = 


0.1540 




1-exp 


8.6( 3)xl0~ 3 


.467( 3) 




2-exp 


1.03( 2)xl0~ 2 


.479 ( 2) 


210/15 


CM 


7.6( 2)xl0~ 3 


.461( 2) 


36/15 


0(a) 


7.6( 2)xl0- 3 


.461( 2) 


36/15 



1.6 



1.4 



1.2 




25 







= 0.1530 




1-exp 


1.05( 3)xl0" 


2 .509( 2) 




2-exp 


1.25( 2)xl0" 


2 .5193(13) 


240/15 


CM 


9.3( 2)xl0~ 3 


.5021(15) 


40/15 


0(a) 


9.3( 2)xl0~ 3 


.5022(15) 


40/15 




K - 


= 0.1455 




1-exp 


3.00( 4)xl0" 


2 .8171(10) 




2-exp 


3.44( 4)xl0" 


2 .8255(10) 


170/15 


CM 


2.47( 5)xl0- 


2 .8073(11) 


44/15 


0(a) 


2.51( 4)xl0" 


2 .8080(11) 


39/15 




K - 


= 0.1420 




1-exp 


4.13( 5)xl0~ 


2 .9576( 9) 




2-exp 


4.70( 5)xl0- 


2 .9655(10) 


136/15 


CM 


3.17( 8)xl0" 


2 .9450(12) 


49/15 


0(a) 


3.30( 7)xl0" 


2 .9466(11) 


41/15 
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= 0.1385 




1-exp 


5.35( 7)xl0" 


2 1.0958( 9) 




2-exp 


6.04( 7)xl0" 


2 1.1033(10) 


115/15 


CM 


3.68(15) x 10" 


" 2 1.0792(14) 


58/15 


0(a) 


4.05(10) x 10" 


- 2 1.0826(11) 


45/15 



Table 2 

Values for x 2 , Z and M for the vector channels 
for a single exponential fit (used as a control), 
a 2-exponential fit, the "QCD Continuum Model" 
(CM) fit to Eq.(0), and the "lattice-QCD Continuum 
Model" (0(a)) fit using Eq.(|). All fitting ranges 
were t = 2 — 20 with the exception of the single ex- 
ponential fit which was t = 12 — 20. 



Figure 1. Plot of the effective mass for the k = 
0.1385 case. The result of a 2-exponential fit to- 
gether with the QCD Continuum Model fits of Eq.(@) 
& Eq.(^) are shown. (The latter two fits are indistin- 
guishable on the plot by eye.) 



cases only: the massless limit ||; and the x = 
limit Jo]. 

In j|] , S q (x) w (for Wilson Quarks) was derived 
as a continuum expansion (i.e. as a Taylor ex- 
pansion in the lattice spacing, a). To obtain this 
expression the authors of M used a special asymp- 
totic series they developed for the modified Bessel 
functions which appear in the standard form of 
S q (x) LAT . Their result, expressed to 0(a) is, 

4m 2 / ram\ r ,\ 3 

I 2 V 



{^fS q {x) w 



0* 2 ) 

2ram Ki 



2ram 



2 . 
m(x 2 ) 
—K\ m(x 2 



(x 2 ) 

+ Am 2 (l + ram)^-^-K 2 

x 2 

+ 0(a 2 ), 



m(x ) ' 



(8) 



where K\ and K2 are modified Bessel functions 
of the second kind. 

In 0, S q (x) w was used instead of Eq.(0) in 
Eq.(|j|) thereby obtaining an expression for G%(t) 
for Wilson fermions correct to 0(a). m lists re- 
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suits for all of the commonly used mesonic chan- 
nels. Here we list only the result for the vector 
channel: 

G Vi {t) OPE,W = _±_ [1 + ra(mi+m2)] 

3 

+ 7T^t[ to i + m 2 - 2raim 2 

+2ra{mim\ + m^)], (9) 
where r is the Wilson parameter. Note that, as 

, j ^,OPE,W ^OPE,cont . ,ru \ 

expected, G 2 ' = G 2 +u{a). 

Table |^ shows the results for the lattice-QCD 
Continuum Model fit, Eq.(||). In this table the 
X 2 values are shown together with the estimates 
of the Z and M values. The lattice-QCD Con- 
tinuum Model fits are very similar to the QCD 
Continuum Model fits from Eq. ((?]), but they be- 
come significantly better (i.e. have a smaller \ 2 ) 
as the quark mass increases. This is to be ex- 
pected since the lattice artefacts ~ ma. (Note 
that in Figure [t] where the effective mass for the 
vector channel at K± t 2 — 0.1385 is plotted, the re- 
sults from the fit to Eq.(||) is indistinguishable by 
eye to that from Eq.(^|).) 

4. Maximum Entropy Method and the 
Nambu Jona-Lasinio Model 

In the above two sections, an Ansatz for the 
spectral function p(s) was derived using an ana- 
lytic expression for the quark propagator, S q (x), 
valid at small x. This Ansatz was then used as a 
fitting function for Monte Carlo data for 2-point 
hadronic functions, G 2 (t). It is interesting to pon- 
der if p(s) itself can be directly extracted from 
Monte Carlo G 2 (t) data. 

At first sight this aim appears impossible to 
achieve due to the intrinsic ambiguity in the in- 
version. G 2 (t) is known at N t w O(10) values of 
t whereas, in order to reveal spectral structure, 
p(s) must be calculated at N s <; 0(100) s val- 
ues. This means that there are an infinite number 
of solutions for p(s) which perfectly fit the data 
G 2 (t). 

In H the Maximum Entropy Method (which is 
well known in almost every numerical field out- 
side lattice gauge theory!) was used to overcome 
this seemingly insurmountable problem. Bayes' 



theorem can be used to express the probability of 
p(s) given the data G 2 (t), P{p\G-i), as 

P(p\G 2 ) <xP(G 2 \p) x P(p). (10) 

Here, P(G 2 \p) can be written as the usual expo- 
nential of the standard \ 2 function of the Maxi- 
mum Likelyhood Method. The probability, P(p) 
can be expressed in terms of the entropy function, 
S. The best solution for p is obtained by "just" 
maximising P(p\G 2 ) in Eq. (|ici|) . 

Since the work of || there have been several 
applications of the MEM method to lattice simu- 
lations. In this paper, we concentrate on applying 
the MEM to the Nambu-Jona-Lasinio model in 
(2+l)D. This work is discussed in more detail in 
H . For details of applying MEM to lattice simu- 
lations see ||]. 

The continuum (Euclidean) NJL Lagrangian is: 

Nt 

The fields and are four-component spinors 
and the index i runs over Nf fermion flavours. 
The model has an interacting continuum limit 
and has been used to model the strong interac- 
tion. It exists in two phases defined by a chiral 
order parameter. The lattice simulation used the 
staggered fermion formulation. 

Figure |^ shows the spectral function obtained 
from MEM for the ir, fermion (f) and massive PS 
meson in the broken phase with Nf — 4 flavours 
and coupling of j3 = 1/g 2 = 0.55. This plot 
represents approximately 30,000 configurations of 
data. It shows that MEM is easily able to isolate 
the ground states in these channels. Furthermore, 
estimates of the binding energy of the PS state 
can be obtained using the MEM spectral func- 
tions which are in agreement with expectations 
from large Nf expansions ||. 

In Figure [|, the ir mass squared is plotted 
against the quark mass showing that the PCAC 
scaling relationship holds. 

Further results for both phases are detailed in 
H and future work where MEM is applied to the 
UKQCD dynamical data set is detailed in llQ] . 
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Figure 2. Broken phase spectral functions for the Figure 3. PCAC scaling relation for (3 = 0.55 NJL 

NJL model. Peaks from left to right are for the tv data sets. 

(m = 0.005, 0.01, 0.02), f and PS (both m = 0.01). 



5. Conclusions 

This paper summarises the results from two 
interesting "unconventional" approaches to the 
study of 2-point hadronic correlation functions in 
lattice Monte Carlo simulations, G-2.it) . The first 
approach aims at obtaining an accurate analytic 
expression for the short distance behaviour of 
C?2(i) via the QCD Sum Rules Continuum Model. 
In particular, the recently derived expression for 
the Wilson quark propagator Q was used. Fur- 
ther work in this area will include the derivation 
of the analytic expression for the clover quark 
propagator and, from that, the calculation of the 
analogous expression for at short distances 

for this action. 

The second approach uses the Maximum En- 
tropy Method to calculate the spectral function, 
p(s), given the 2-point function, G*2(i). This tech- 
nique is being used by a growing number of re- 
searchers in this context. In this work we used 
MEM to study the Nambu-Jona-Lasinio model 
and found it to work admirably. 
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